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1. Introduction 

The twisted mass formulation of Lattice QCD [1, 2] is being studied extensively with Nf = 2 
dynamical flavours by the European Twisted Mass (ETM) collaboration [3-7]. In this formulation 
of QCD, the Wilson term is chirally rotated within an isospin doublet. To include a dynamical 
strange quark in a unitary setup, we add, in addition to the strange quark a charm quark in a heavier 
and mass-split doublet as discussed in [8-10]. We will briefly describe our action in section 

2, recapitulate our procedure for tuning to maximal twist and focus on the tuning of the heavy 
doublet. We give an overview of the runs we have carried out and section 3 gives first results for 
some light-quark sector observables. 

2. Lattice setup 

In the gauge sector we use the Iwasaki gauge action [1 1]. With this gauge action we observe a 
smooth dependence of (possible) phase sensitive quantities on the hopping parameter K around its 
critical value ff cr j t . The fermionic action for the light doublet is given by: 

Si = a 41 £{Xiix)[Dw[U]+m j + iiiiY5T3]Xi(x)}, (2.1) 

X 

using the same notation as used in [10]. In the heavy sector, the action becomes: 

S h = a 4 Y,{%h{x) [D W [U] +mo,h + /MaTSTi + Ms^] Xh(x)} . (2.2) 

X 

At maximal twist, physical observables are automatically 6(a) improved without the need to de- 
termine any action or operator specific improvement coefficients. The gauge configurations are 
generated with a Polynomial Hybrid Monte Carlo (PHMC) updating algorithm [12-14]. 

2.1 Tuning action parameters 

Tuning to maximal twist requires to set m^j and mqj, equal to some proper estimate of the 
critical mass m CI n = m cr i t (j3) [8]. Here we set m j = m^/, = I /(2k) —4. As has been shown in 
[9], this is consistent with 0(a) improvement defined by the maximal twist condition ampcAC.i = 
(see also ref. [10]). The numerical precision at which the condition nipcACj = is fulfilled in order 
to avoid residual large &(a) effects when the pion mass is decreased is, for the present range of 
lattice spacings, |e//i/| < 0.1, where e is the deviation of rapcAC./ from zero [4, 15]. As explained 
in [10], tuning to K^t was performed independently for each \i\ value. From table 1 we observe 
that the estimate of K^t depends weakly on The heavy doublet mass parameters \i a and Us 
should be adjusted in order to reproduce the values of the renormalized s and c quark masses. The 
latter are related to jj, a and jUg via [8] : 

1 Z P 

(m SiC ) R = —(n a =f—n s ), (2.3) 

where the — sign corresponds to the strange and the + sign to the charm. In practice we fix the 
values \i a and \i§ by requiring the resulting K- and D-meson masses to match experimental results. 
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Label 


P 


K 


ajii 


aji a 




L/a 


T/a 


m n L 


I £ /m/| 


C\ 


1.90 


0.1632700 


0.0040 


0.150 


0.190 


20 


48 


3.0 


0.14(14) 


c 2 


1.90 


0.1632700 


0.0040 


0.150 


0.190 


24 


48 


3.5 


0.07(14) 


Ai 


1.90 


0.1632650 


0.0060 


0.150 


0.190 


24 


48 


4.1 


0.03(3) 


A 2 


1.90 


0.1632600 


0.0080 


0.150 


0.190 


24 


48 


4.8 


0.02(2) 


AjDi 


1.90 


0.1632550 


0.0100 


0.150 


0.190 


24 


48 


5.3 


0.02(2) 


A 4 


1.90 


0.1632720 


0.0030 


0.150 


0.190 


32 


64 


4.0 


0.08(7) 


A5C3 


1.90 


0.1632700 


0.0040 


0.150 


0.190 


32 


64 


4.5 


0.04(5) 


A 6 


1.90 


0.1632670 


0.0050 


0.150 


0.190 


32 


64 


5.0 


0.05(2) 


D 2 


1.90 


0.1632550 


0.0100 


0.150 


0.197 


24 


48 


5.3 


0.35(1) 


Bi 


1.95 


0.1612400 


0.0025 


0.135 


0.170 


32 


64 


3.4 


0.06(6) 


B 2 


1.95 


0.1612400 


0.0035 


0.135 


0.170 


32 


64 


4.0 


0.02(2) 


B 3 


1.95 


0.1612360 


0.0055 


0.135 


0.170 


32 


64 


5.0 


0.08(1) 


B 4 


1.95 


0.1612320 


0.0075 


0.135 


0.170 


32 


64 


5.8 


0.05(1) 


B 5 


1.95 


0.1612312 


0.0085 


0.135 


0.170 


24 


48 


4.6 


0.01(2) 



Table 1: Input parameters, m K L and | e / ju./ 1 for all ensembles used in this paper. Every ensemble has ~ 5000 
thermalized trajectories of length z = 1. We have two main ensemble sets: A and B, at J3 = 1.90 and 
j3 = 1.95 respectively. Ensembles labeled C are used to check finite size effects. Ensembles labeled D are 
used to check/tune the strange and charm quark masses. 

2.2 Determination of heavy-light meson masses 

Since the twisted mass lattice Dirac operator of the non-degenerate heavy quark doublet (cf. 
(2.2)) contains a parity odd and flavour non-diagonal Wilson term, parity as well as flavour are 
not anymore quantum numbers of the theory. In contrast to parity and flavour conserving lattice 
formulations, it is not possible to compute correlation functions restricted to a single parity and 
flavour sector in this setup. While the /T-meson will remain the lightest state and therefore relatively 
easy to extract, for a theoretically clean determination of the D-meson mass one has to consider the 
four sectors labeled by parity = ± and flavour = sjc at the same time. And since besides the 
^-meson there are a number of K + n x 71 states and possibly also "positive parity K states" below 
the D-meson, this renders the D a highly excited state. At currently available statistics it seems 
extremely difficult to extract such a high lying state. 

As such, we resort to a different strategy in order to extract this mass. We attempt to determine 
the mass of the D-meson without computing the full low-lying spectrum, e.g. we do not determine 
all low lying states below the D. To this end we apply smearing techniques (cf. [16], where the 
same setup was used) to construct highly optimized trial states with large overlap to the K- and D- 
meson, and make certain assumptions about these trial states, which will be motivated and detailed 
in an upcoming publication. To extract the D-meson mass, we applied three different methods: (1) 
solving a generalized eigenvalue problem, (2) performing a multi-exponential fit and (3) rotating 
the twisted basis correlators back to the physical basis (in order to do this we need to compute the 
light and heavy twist angle and a ratio of the appropriate renormalization constants). The values of 
the D-meson mass extracted from these three different methods are consistent with each other. 
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2.3 Status 

The left panel of figure 1 shows the tuning of the strange quark mass by showing the difference, 
scaled with the chirally extrapolated value of ro/a between twice the T^-meson mass squared and 
the pion mass squared. Set A at /3 = 1.90 a\ig = 0.190 (green points) appears to overshoot the 
physical point (the black cross on the left), while set B (red points) extrapolates better. To improve 
the tuning of the strange quark mass for set A, we are currently applying a reweighting procedure 
as described in [10] in the parameters apLg and K. The blue point with a different heavy sector 
splitting apLg = 0.197 is a run to check this procedure. Though this run is not tuned to maximal 
twist yet, the i^-meson mass appears to be much closer to its physical value. The right panel of 
figure 1 shows the mass of the D-meson (obtained in this case by method (3)) as a function of the 
pion mass squared for various simulation points as well as the experimental value from the Particle 
Data Group [17] . The plot demonstrates that we have tuned the charm (sea) quark mass in our 
simulations to a physically realistic value. As a final check, we also use an estimate of Zp/Zs to 
verify that m c ~ \Qm s . 
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Figure 1: r\(2ni\ — mfy and r^mu as functions of (rom K ) 2 , showing the status of the tuning of the sttange 
and charm quark mass respectively. The experimental value from PDG is added as the black cross (/n = 
0.45(3) fm was used). Red points label the j5 = 1.95 runs, green points label the j3 = 1.90 runs, where the 
single blue point corresponds to j3 = 1.90 with a different heavy sector splitting ajig. Circles denote runs 
with L/a = 24, triangles indicate a volume with L/a = 32. 



3. Results 

As a first check of our data, we have compared it to the extensively analysed data set that exists 
for our Nf = 2 data. To compare the two sets, we plot dimensionless physical ratios in figure 2. The 
figure shows no evidence of disagreement between all our results, suggesting small discretisation 
effects and small effects of dynamical s- and c-quarks for these observables. 
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Figure 2: (m % j f K ) 2 vs [m K /m N ) 2 (left) and r / s vs (r Q tn K ) 2 (right) for both A/y = 2+ 1 + 1 data (with 
j3 = 1.90, 1.95) and Nf — 2 data (with /3 = 3.9,4.05,4.2 and using a different gauge action). In the left 
plot, niN is the nucleon mass, and the physical point is included as the black cross. In both plots finite size 
corrections are not applied, and in the right plot the chirally extrapolated values for ro/a were used. 



3.1 Light meson chiral perturbation theory fits 

In order to extract the lattice spacing and light quark mass from our data-sets, we perform a 
next to leading order SU (2) chiral perturbation theory fit of the m % and f n data. We use continuum 
formulae and correct for finite size effects either without any new low energy constants a la Gasser 
and Leutwyler [18], or with l\ and T% added in, as described in [19]. The results are listed in 
table 2. We have performed these fits for ensemble sets A and B separately, and also combined 
them in a single fit. In table 2, we include a systematic error, estimated at 2 — 5%, coming from 
the dispersion of the values of the fitted parameters between NLO and NNLO. Note that since the 
quark mass enters the #PT expression, in order to combine the two sets at different lattice spacings, 
we need to know the renormalization factor of the quark mass = 1 /Zp, a computation which is 
not yet complete. Assuming that Zp is effectively a function of /3 in the range of parameters we are 
considering, we can fit the ratio of those Zp-values and lattice spacings and extract lattice spacings 
from the combined fit. In every fit we use as inputs the physical f n and m % , and extract fo, h and 
I4. A complete analysis (analogous to [20]) of the systematic effects is in progress. 



set 


pts 


/o(MeV) 


h 


u 


«j3=1.90( fm ) 


«|3=1.95( fm ) 


A&B 


11 


121(4) 


3.5(2) 


4.7(2) 


0.086(6) 


0.078(6) 


A 


6 


121(4) 


3.4(2) 


4.8(2) 


0.086(7) 




B 


5 


121(4) 


3.7(2) 


4.7(2) 




0.078(7) 



Table 2: Results from the NLO SU (2) jPT fits for combined, only set A and only set B respectively. Errors 
are dominated by a systematic error of 2 — 5% due to performing an NLO fit. The column "pts" refers to the 
number of ensembles used in that fit. 
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3.2 Chiral extrapolation of the nucleon mass 

In this section, we present preliminary results for the light quark mass dependence of the 
nucleon mass. We consider the one-loop result from heavy baryon chiral perturbation theory 
(HB^PT) 

m N = m%- Acxm\ - — A m 3 „ (3.1) 

and fix the scale and light quark mass to the point where the ratio m^/m n attains its physical value. 
We fix f n and gA to their physical values (130.7 MeV and 1.27 respectively) as has also been 
performed in [7]. Using this procedure, we find a lattice spacing of 0.089(2) fm and 0.077(3) 
fm for j3 = 1.90 and 1.95 respectively. The % 2 /(d.o.f.) of these fits is not very good, and fitting 
a linear extrapolation appears to be consistent with the data. This is not unique to our data, and 
has been observed by various collaborations. We therefore perform the linear fit here as well, 
and absorb the difference between the two extrapolations in the systematic error. A more detailed 
analysis of the chiral extrapolation of the nucleon mass will be presented in an upcoming study. 
The lattice spacings that we obtain from the chiral extrapolation of the nucleon mass are 0.089(9) 
fm and 0.077(4) fm for set A at j8 = 1.90 and set B at /3 = 1.95 respectively. 

3.3 ro / a extrapolation 

Since ro/a is very sensitive to K in the vicinity of fQ r jt, the fact that we now tune to maxi- 
mal twist at every value of /I/, might, with respect to what was done for the Nf = 2 case, in part 
provide an explanation for the observed change of slope in the mass dependence of ro/a between 
Nf = 2 and Nf ■ = 2 + 1 + 1. Note however that these differences tend to diminish when increas- 
ing the value in j3 in the Nf = 2 + 1 + 1 case. We extrapolate ro/a using a simple quadratic fit 
ro/a = c\ +C2« 2 i u / 2 , where c\ is the value of ro/a in the chiral limit. We perform both a polyno- 
mial fit ro/a = c\ +C2afii +c^a 2 pL } 1 and a linear fit ro/a = c\ -\-cia\Li to help estimate systematic 
errors. We find that based on the ^ 2 /d.o.f. the quadratic fit is for both values of j3 the best fit. The 
polynomial fit gives nearly identical results for c\, while c\ from the linear fit is 1 to 3a higher. 
Using the lattice spacings from the combined light meson chiral perturbation theory fit, we extract 
two predictions for ro, which seem to agree well at ro = 0.45(3) fm. 



p 


c\ (quadratic fit) 


a(fm) 


r (fm) 


1.90 
1.95 


5.24(2) 
5.71(4) 


0.086(6) 
0.078(6) 


0.45(3) 
0.45(3) 



Table 3: ro determination for both ensembles separately, cj(qua) is the value of a quadratic ro/a extrap- 
olation in the chiral limit with the statistical error in brackets. The lattice spacings a are taken from the 
combined light meson chiral perturbation theory fit. The obtained values for ro from the two ensembles 
seem to agree well with each other. 

4. Conclusions 

We have presented first results from runs performed with Nf = 2+ 1 + 1 flavours of dynamical 
twisted mass fermions. No evidence of disagreement between these results and those with Nf = 2 
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twisted mass fermions is shown through dimensionless ratio plots (of m n , f K , and ro/a), sug- 
gesting small discretisation effects and small effects of dynamical s- and c-quarks for these observ- 
ables. We have extracted the lattice spacings of our two ensemble sets using two different methods, 
which agree within errors with each other. We have measured ro on both ensembles and found 
consistent results. We are in the process of performing a more detailed combined analysis in order 
to improve our understanding of the systematic errors. 
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